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Abstract 

In this work we report the results concerning calculations for quantum-mechanical rotational 
transitions in molecular hydrogen, H2, induced by an ultra-cold ground state anti- hydrogen atom 
His. The calculations are accomplished using a non-reactive close-coupling quantum- mechanical 
approach. The H2 molecule is treated as a rigid rotor. The total elastic scattering cross section 
(7e;(e) at energy e, state-resolved rotational transition cross sections o-jji{e) between states j and j' 
and corresponding thermal rate coefficients kjj'{T) are computed in the temperature range 0.004 
K ^ r < 4 K. Satisfactory agreement with other calculations (variational) has been obtained for 
o-ei(e). 
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I. INTRODUCTION 



Interaction and collisional properties between matter and antimatter is of fundamental 
importance in physics l|, l2|. The anti- hydrogen atom H, which is a bound state of an anti- 
proton p~ and a positron e"*", is the simplest representative of an antimatter atom. This is a 
two-particle system, which, however, may possess very different interactional and dynamical 
properties compared to its matter counterpart the H atom 4|. 

By now much effort has been exerted in various experiments to build and store H at cold 
and ultra-cold temperatures [sl, I^, Q, I^, [o] . New experiments are planned or in progress to test 
the fundamental laws and theories of physics involving antiparticles and antimatter in general 
[2I. For example, it follows from the CPT symmetry of quantum electrodynamics that a 
charged particle and its anti-particle counterpart should have equal and opposite charges and 
equal masses, lifetimes and gyromagnetic ratios. The CPT symmetry predicts that hydrogen 
and anti-hydrogen atoms should have identical spectra. In this way experimentalists plan 
to test whether in fact H and H have such properties. Specifically, a starting point would 
be to compare the frequency of the ls-2s two-photon transition in H and H. Also, one of 
the important practical applications of antihydrogen has been mentioned in [l^, where the 
authors considered controlled H propulsion for NASA's future plans in very deep space. 
Researchers at CERN jl| and from other groups 8] are interested to trap and study H at 
low temperatures, e.g., T < 1 K, when the H atom will be almost in its rest frame. The 
study of Lamb shift and response of antihydrogen to gravity at ultra-low energies should 
allow them to test more precisely the predictions of two fundamental theories of modern 
physics: quantum field theory and Einstein's general theory of relativity 

It has been pointed out that the main cause of loss of H atoms confined in a magnetic 
gradient trap is due to H-I-II2 and H+He collisions. Therefore, the H-I-H2 scattering cross- 
sections and corresponding rotational- vibrational thermal rate coefficients, in the case of H2, 
would be very helpful to gain a practical understanding of the slowing down and trapping 
of H. Hence the investigation of the possibility of cooling of H atoms by colliding them with 
colder H2 is of significant practical interest 12|. (Similar collision between trapped fermionic 
atoms with cold bosonic atoms has been fundamental in cooling the fermionic atoms and thus 
leading them to quantum degeneracy Such investigation of H interaction with H and 

H2 can reveal the survival conditions of H in collisions with H and, even more importantly. 
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with H2 Q. 

Further, coohng occurs by energy transfer in elastic colhsions of H with H2. However, 
during the colhsion the rearrangement process may lead to the formation of protonium {pp~) 
and positronium {e~^e~) exotic atoms and the destruction of H atoms. They are formed as 
matter-antimatter bound states, which then annihilate. (There have also been many studies 
of scattering of positronium atoms 15|], the lightest matter-anti- matter atom). Thus one 
can conclude, that the effectiveness of cooling H is determined from a comparison of the 
cross sections for direct scattering and rearrangement. 



By now a series of theoretical works have been published, in w 



interaction between H and H, He, H2 have been investigated 



16, 



17,hLl 



lich the properties of 
□ 

1191 ] . Some theoretical 



studies have been carried out for the h 
mechanical methods l20 



21 
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H system at thermal energies using quantum- 



251 . Also, discussions on the importance and 



applications for this system, especially, in connection with Bose-Einstein condensation 26| , 
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21| , and its static and dynamic properties 27|| , can be found in the 



ultra-cold collisions 
literature. 

In this work we present results for the collision of an ultra-cold H atom with H2, where 
H2 is treated as a rigid rotor with a fixed distance between hydrogen atoms. The elastic, 
rotational state-resolved scattering cross sections for the H-H2 scattering and their corre- 
sponding thermal rate coefficients are calculated using a non-reactive quantum-mechanical 
close-coupling approach. The potential interaction between H and the hydrogen atoms is 



taken from Ref. 



16|- 



In the next section we present the quantum-mechanical formalism used in this work. The 
results and discussion are presented in Sec. HI. Conclusions are given in Sec. IV. 



II. H-H2 SCATTERING FORMULATION 



A. Basic Equations 

In this section we describe the close-coupling quantum-mechanical approach we used to 
calculate the cross sections and collision rates of a hydrogen molecule H2 with an anti- 
hydrogen atom H. Atomic units {e = nie = h = 1) are used in this section, where e and 
are charge and mass of an electron. Three-body Jacobi coordinates {r, i?} for the H-|-H2(j) 
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system used in this work are shown in Fig. 1. The two H atoms are labeled 2 and 3 and the 
H atom is labeled 1, O is the center of mass of the H2 molecule, is the polar angle between 
vector r connecting the two H atoms in H2 (labeled 2 and 3) and vector R connecting the 
center of mass of the H2 molecule to the H atom (labeled 1). Next, j and L are angular 
momenta corresponding to the vectors r and R, respectively. The quantities X21 and X31 are 
the distances between the H atom labeled 1 and the H atoms labeled 2 and 3, respectively. 

The Schrodinger equation for an a + 6c collision in the center of mass frame, where a (H) 
is an atom and be (H2) is a linear rigid rotor, is 
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^ + l^ + i^W^)-^] = 0- (1) 

where is the relative momentum between a and be, Mr is the reduced mass of the atom- 
molecule (rigid rotor in this model) system a + be: Mr = rriaimb + rric) / {ma + nib + nic), /i 
is the reduced mass of the target: /i = marrih/ {ma + ?Tife), f is the angle of orientation of the 
rotor ab, V{r, R) is the potential energy surface (PES) for the three-atom system abe, and 
E is the total energy of the system. The eigenf unctions of the operator L^? in Eq. ([T]) are 
the spherical harmonics Yjm{T) 



To solve Eq. ([Tj), the following expansion is used [30| 

R 



^{r,R)= Y: ^^^0/f(r,^), (2) 



JMjL 

where channel expansion functions are 



(/./f (f, R)=J2 C^ZLmJjm.Ar)YLm,{R), (3) 



lJM, 

mim2 

here J = j + L is the total angular momentum of the system abe, and M is its projection 
onto the space fixed z axis, mi and m2 are projections of j and L respectively, C'-^^j^^^ are 
the Clebsch-Gordan coefficients, and f/'s are the appropriate radial functions. 

Substitution of Eq. ([2]) into Eq. ([T]) provides a set of coupled second order differential 
equations for the unknown radial functions Ujl^{R) 

^ - ^^^^ + ^Jl) U^^'{R) = 2MrY. I < <p^t\f,R)\V{r,R)\^fl{f,R) > 

xU}l^,{R)dfdR. (4) 

To solve the coupled radial equations (jlj), we apply the hybrid modified log-derivative- Airy 
propagator in the general purpose scattering program MOLSCAT |3l|. Additionally, we 



tested other propagator schemes included in MOLSCAT. Our calculations reveal that other 
propagators can also produce quite stable results. 

The log-derivative matrix is propagated to large intermolecular distances R, since all ex- 
perimentally observable quantum information about the collision is contained in the asymp- 
totic behavior of functions U^l^{R oo). The numerical results are matched to the known 
asymptotic behavior of Uj^'^{R) relating to the the physical scattering S'-matrix 32 1 



where = [2Mr{E - E^)]^/'^ is the channel wave-number of channel a = (jL), E^ is 
rotational channel energy and E is the total energy in the ahc system. This method was 
used for each partial wave until a converged cross section was obtained. It was verified 
that the results have converged with respect to the number of partial waves as well as the 
matching radius, Rmaxi for all channels included in our calculations. 

Cross sections for rotational excitation and relaxation phenomena can be obtained directly 
from the S'-matrix. In particular, the cross sections for excitation from j — > j' summed over 



the final m' and averaged over the initial m are given by 30 1 

TT 

(2j + l)kl 



= Y.^'iJ +mn'hL' - S\fL'-jL-E)\\ (6) 



JLV 

The kinetic energy is e = E — Bej{j + 1), where Be is the rotation constant of the rigid rotor 
be, i.e. the hydrogen molecule. 

The relationship between the rate coefficient kj^j>{T) and the corresponding cross section 



aj^j'{Ekin) can be obtained through the following weighted average [33 1 



k^-AT) = ^TT^ / a,^Ae)e-^'''^ede, (7) 



-kMr {kBTf 

where e = E — Ej is pre-collisional translational energy at temperature T, ks is Boltzman 
constant, and is the minimum value of the kinetic energy needed to make Ej levels 
accessible. 
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B. H— H2 Interaction Potential 



In the following section, we will present our results for rotational quantum transitions in 
collision between H2 and an anti-hydrogen atom H, that is 

H2(j) + H^H + H2(/). (8) 

Here H2 is treated as a vibrationally averaged rigid monomer rotor. The bond length was 
fixed at 1.449 a.u. or 0.7668 A. The rotation constant of the H2 molecule has been taken as 
Bp = 60.8cm~^. The Ho rigid rotor model has been already applied in different publications 



3G 
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401 ]. For the considered range of kinetic energies the model 
can be quite justified in this special case when only pure rotational quantum transitions 
at low collisional energies are of interest as in H2(j)+H, and when the energy gap between 
rotational and vibrational energies is much larger than kinetic energy of collision. In such a 
model the quantum mechanical approach is rather simplified. 

Next we consider an important physical parameter in atomic and molecular collisions, 
e. g., the PES between the atoms. There is no global potential energy surface available 
for the three-atom H-H2 system. However in Ref. jl6], the author calculated the values of 
interaction energy between H and H, i. e., the H-H energy curve using the Rayleigh-Ritz 
variational method. Further, the microHartree accuracy of Born-Oppenheimer energies of 
the system has been achieved in that work. 

To construct the H2-H interaction potential we take the H-H energy data from Ref. 
and make a cubic spline interpolation through all 46 points taken from Table 1 of that paper. 
These data have been tabulated from Rmin = 0.744 a.u. to Rmax = 20.0 a.u. inter-atomic 
distances. Because in the current work we use the rigid rotor model for H2, we do not need 
the interaction energy between hydrogen atoms in H2. The interaction potential between a 
hydrogen molecule and H is taken by sandwiching the two H-H potential energy curves: 

V{r, R) = V{r, R, 6) = V^\j^{x2i) + V^lj^{x3^), (9) 

where distances between atoms are written as follow (cf. Fig. 1): 

X21 = Vr2/4 + i?2 + ^i?cose and X31 = v^r2/4 + R^ + ri? cos(7r - 6). (10) 

The functions V^^^{y) with k = 2(3) are represented as cubic spline interpolation func- 



tions for any value of y = X21 or y = X31 as follows: 

Vi'_M = ^H-h(^^) + B^(y - + C^{y - X,) + D,{y- X,), (11) 

where Bi{y — Xi), Ci{y — Xi), Di{y — Xi) perform the spline interpolation and where ^ 
y ^ Xi+i, in each sub-interval [Xi, Xj+i], i = 1, 2, 3, (n — 1), n = 46. The coordinates Xi 
and corresponding values of the H-H potential energy data have been taken from Table 1 of 
Ref. [161 . 



The calculated PES is shown in Fig. 2. It is clear that the potential has a singular 
value when the distance between H and H2 is equal to zero. Additionally, the H-H2 PES 
which we obtain from Eq. is shown in Fig. 3. Specifically, this potential was used in our 
calculations of H + H2 collisions. Again, as seen in Fig. 2 the potential energy curve between 
H-H has a Coulomb type singularity at small distances. In our calculations we needed to 
make additional test runs to achieve convergence in our results. In the next section we 
will briefly demonstrate the numerical convergence of the results when calculating total 
elastic scattering cross sections. These results depend on various numerical and quantum- 
mechanical scattering parameters. 



III. RESULT AND DISCUSSION 



A. Convergence test 

Numerous test calculations have been undertaken to insure the convergence of the results 
with respect to all parameters that enter in the propagation of the Schrodinger equation. 
These include the atomic-molecular distance R, the total angular momentum J, the number 
of total rotational levels to be included in the close-coupling expansion and others, see Fig. 1. 
Particular attention has been given to the total number of numerical steps in the propagation 
over the distance R of the Schrodinger equation (jlj). Speciflcally, the parameter R ranges 
from 0.75 a.u. to 20.0 a.u. We used up to 50000 propagation points. We also applied and 
tested different mathematical propagation schemes included in MOLSCAT. 

The rotational energy levels of para-B.2{j) and the corresponding angular momenta j 
are shown in Table I. The goal of this work is to get new results for H + para-'H.2 thermal 
rate coefficients kj^j>{T) at ultra-low temperatures: specifically 0.004 K< T < 4 K. The 
corresponding cross sections have been calculated for collision energies varying from ~ 0.0001 
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cm~^ to ~ 100 cm~^. These energy values are very small. However, despite this fact, to 
reach convergence of the results we needed to include in expansion ([2]) a significant number 
of rotational levels of the H2 molecule, specifically up to jmax ~ 60. Below in Table II we 
present these results. 

Also, we found that, at lower energies, for the numerical solution of Eq. (4) a much larger 
number of propagation (integration) points are needed than at higher energies. Specifically, 
at higher energies we need 500 propagation points, but for lower energies 50,000 points are 
needed to achieve comparable precision. Convergence has been achieved for elastic scattering 
cross sections for various scattering parameters. Below in Table III we present these results. 
Then we used this data in our calculation for rotational energy transfer, elastic scattering 
cross section and thermal rate coefficient. 

Now in Table II we present results for the total elastic cross sections for two coUisional 
energies: 0.1 cm~^ and 0.01 cm~^ The cross sections are shown for a number of different 
maximum values of the rotational angular momentum j = j^ax in the H2 molecule included 
in expansion ([2]). This is also the JMAX parameter in MOLSCAT 3l|- Other scattering 
parameters have also been treated correctly in the calculation. One can see that JMAX 



should be at least 56. The other scattering parameter in Table II is MXSYM 
the number of terms in the potential expansion over angular functions 



29 



3l|. It reflects 



301 . It is therefore 



evident that in this calculation we need to keep at least 24 terms in the expansion. 

In Table III we also present results for total elastic scattering cross sections for few more 
selected energies. However in this table the convergence has been reached by increasing the 
total angular momentum J and by increasing the total number of the propagation steps in 
the propagation over the coordinate R of the Schodinger equation (j4]). As expected, for lower 
energy collision we needed smaller values for the maximum J. For example, for collision 
energy Ecoii = 0.01 a.u. it is enough to have J = 0, however J = 10 should be taken for 
E,ou = 100.0 a.u. 

In regard to the total number of the propagation steps in the solution of Eq. (j4]) one 
can see, that we need to include many more propagation points for low energy calculations. 
All test calculations in Table HI have been done with JMAX = 56 and MXSYM = 24. 
The obtained results concerning the numerical and scattering parameters have been used in 
our calculation for total elastic aei{E) and rotational quantum state transfer aj^j>{E) cross 
sections and corresponding thermal rate coefficients kj^ji{T). 
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B. H + para-'H.2 results 

Now we present computational results for process ([8]), namely, for elastic scattering {j = 
^ j' = 0) and for low quantum number rotational transition between levels with j = 0, 2 
and 4: 2 — > 0, — > 2, 4 — 0, and 4 — 2. From the results of Table III we see that to 
reach numerical convergence, for example, for the elastic scattering cross section, we need 
to include a large number of H2 rotational levels, specifically up to 60. 

The results for the elastic scattering cross sections a^i (e) for H + H2 — H2 + H are shown 
in Fig. 4 tog ether with the corresponding results of variational calculations of Gregory and 



Armour [iJ]. It can be seen, that basically the two sets of cross sections are close to 
each other, although in our calculation we use larger number of collision energy points, 
specifically up to 200. In our calculation a shape resonance is found at energy e ~ 3.5 x 10~^ 
Hartree. As in Ref. [l^ our o"e;(e) tends to reach a constant value at lower energies with 
(Te;(e < 10~® a.u.) = 9.47 x lO^a^. This result allows us to calculate the H+H2 scattering 
length, which is 

a = Vae//(47r) = 27.5 ao- (12) 

The Gregory-Armour scattering length jl^ obtained with a variational method is a = 
19.5 ao. The two results are in reasonable agreement with each other. 

In Fig. 5 (a) and (b) we show the total state-resolved cross sections aj=2->j'=o{'v) vs. 
velocity v and the corresponding thermal rate coefficients A;2^o(T) vs. temperature T for the 
hydrogen molecule rotational relaxation process. It is seen, that when the collision energy in- 
creases the de-excitation cross section decreases. It can be explained in the following way: at 
low relative velocities (kinetic energies) between Il2(j = 2) and H, the rotationally excited H2 
molecule has more time for interaction and consequently, it has higher quantum-mechanical 
probability to release its internal rotational energy to H. The resulting corresponding rate 
coefficients have been calculated for a temperature range from 0.004 K < T < A K and 
are also presented in Fig. 5 below the cross section results. 

Next, in Fig. 6 (a) and (b) we present results for the total state- resolved cross sections 
(To^j'=2(v) vs. velocity v and the corresponding thermal rate coefficients A;o^2(T) vs. temper- 
ature T for the hydrogen molecule rotational excitation process. It is quite understandable, 
as we find from Fig. 6 (a), that when the collision energy (relative velocity v) is increases the 
quantum-mechanical probability and corresponding cross section of the rotational excitation 
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of H2(j) are also increases. Figure 6 (b) depicts the corresponding results for the thermal 
rate coefficient. 

In Figs. 7 (a) and (b) we present results for cross sections and rates for rotational 
relaxation process, as in Figs. 5 (a) and (b), but now connecting the states j = 4 and j' = 2. 
Finally, in Figs. 8 (a) and (b) we present results for cross sections and rates for rotational 
relaxation process but now connecting the states j = 4 and j' = 0. An unexpected result 
has been found in Fig. 7 (a) in the rotational transition cross section (74-^2 (f), i-e., when the 
H2 quantum angular momentum has been changed from j = 4 to j' = 2. One can see, that 
the values of these cross sections at very low collision energies are almost from 5 to 10 times 
larger then other cross sections considered in this work, compare with the results from Figs. 
5, 6 and 8. 



IV. SUMMARY 



A quantum-mechanical study of the state-resolved rotational relaxation and excita- 
tion cross sections and thermal rate coefficients in ultra-cold collisions between hydrogen 
molecules H2 and anti-hydrogen atoms H has been the subject of this work. A model PES 
for H2-H has been constructed by sandwiching two H-H interaction potentials for two differ- 

n - 

ent hydrogen atoms taken from Ref. pj|. This H-H interaction potential is shown in Fig. 2. 
The H2-H PES is presented in Fig. 3. Calculation for total elastic scattering cross section 
and for low quantum rotational transition states have been performed. We considered only 
the following quantum transitions: 2^0, 0^2, 4^2, and 4^0. 

A test of the numerical convergence was undertaken. These results are presented in Tables 
II and HI. Our results reveal that it is necessary to set the rotational angular momentum 
jmax in the H2 molecule to a relatively large number, i.e. in the expansion ([2]) we needed 
io include up to 60 terms. The calculation was performed using the MOLSCAT program 



3l|. Different propagation schemes included in the MOLSCAT program have been used 



and tested. Additionally, the MXSYM potential parameter in that program also needed 
to have a relatively large value to obtain good convergence as can be seen from Table II. 
The numerical convergence has also been tested over the total number of the propagation 
steps over coordinate R in the solution of the Schrodinger equation (jlj). We have found, 
that at low energies we need a much larger number of integration points than at higher 
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energies, cf. Table III. Our results for the H2(j)+H total elastic scattering cross section are 
in reasonable agreement with the corresponding results from Gregory and Armour jl^ . The 
authors of this paper used a different PES, which is still unpublished, and applied a quantum- 
mechanical variational approach. Unfortunately, the rotational transitions in the H2(j)+H 



3- 



collisions have not been calculated in that work [IJl. One of the interesting results of the 
present work is that the cross section of the rotational transitions from H2(j = 4 ^ j = 2) 
at ultra-low energies are approximately 5-10 times larger than other transition state cross 
sections. 

To the best of our knowledge we do not know of any other calculation of the rotational 
transitions in the H2(j)-|-H collision. These results can help to model energy transfer pro- 
cesses in the hydrogen-anti-hydrogen plasma, and perhaps to design new experiments in the 
field of the anti-hydrogen physics. Finally, we believe, that in the future work it should 
be useful to include vibrational degrees of freedom of the H2 molecules, i.e. to carry out 
quantum-mechanical calculations for different rotational-vibrational relaxation processes: 
B.2{v,j) -|- H — > H -|- B.2{v',j'), where v and v' are the vibrational quantum numbers of H2 
before and after the collision, respectively. 
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TABLE I: para-il2 rotational spectrum. 



Level Rotational energy (cm ^) Internal quantum momentum in para-H2(j) 

1 0.00 

2 364.80 2 

3 1216.00 4 

4 2553.60 6 

5 4377.60 8 

6 6688.00 10 

7 9484.80 12 

8 12768.00 14 

9 16537.60 16 

10 20793.60 18 

11 25536.00 20 

12 30764.80 22 

13 36480.00 24 

14 42681.60 26 

15 49369.60 28 

16 56544.00 30 

17 64204.80 32 

18 72352.00 34 

19 80985.60 36 

20 90105.60 38 

21 99712.00 40 

22 109804.80 42 

23 120384.00 44 

24 131449.60 46 

25 143001.60 48 

26 155040.00 50 

27 167564.80 52 

28 180576.00 54 

29 194073.60 56 
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TABLE II: Convergence of the total elastic scattering cross section cjez (10~^^cm^) at different 
energies E (cm~^) in II+II2 — > H2+II with respect to the maximum value of the rotational angular 
momentum j = jmax in H2(i) included in the expansion Q (parameter JMAX in MOLSCAT). Con- 
vergence with the number of terms in the potential expansion (parameter MXSYM in MOLSCAT) 
is also shown. Numbers in parentheses are powers of 10. 



E (cm-i) 










(Je^X 10^6 


(cm^) 












JMAX 








MXSYM 






30 


40 


50 


56 


60 


12 


20 24 


26 


0.1 
0.01 


61.0 
55.4 


5.25(2) 
1.06(3) 


1.59(3) 
6.59(3) 


1.59(3) 
6.62(3) 


1.59(3) 
6.62(3) 


1.97(3) 
1.12(4) 


1.61(3) 1.59(3) 
6.75(3) 6.56(3) 


1.59(3) 
6.54(3) 
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TABLE III: Convergence for the total elastic scattering cross section aei (10~^^cm^) at different 
collision energies E in (cm~^) in H+H2 H2+H with respect to the maximum value of the total 
angular momentum J of the 3-atomic system: parameter JTOT in MOLSCAT. Convergence on 
the number of numerical space steps in propagation over distance R of the Schrodinger equation 
(parameter STEPS in MOLSCAT) is also shown. Numbers in parentheses are powers of 10. 



E (cm^i) 








< 10^6 (cm^) 












JTOT 






STEPS 











2 




500 


1000 


10000 


50000 


0.1 


1.59(3) 


1.59(3) 




3.64(2) 


2.67 


1.60(3) 


1.59(3) 


0.01 


6.62(3) 


6.62(3) 




3.62(2) 


1.75(-1) 


6.63(3) 


6.62(3) 




4 


6 


8 


500 


1000 


5000 


7000 


10.0 


5.85(1) 


5.96(1) 


5.96(1) 


6.23(1) 


6.00(1) 


5.96(1) 


5.96(1) 


1.0 


1.69(2) 


1.69(2) 




1.77(2) 


1.70(2) 


1.70(2) 


1.70(2) 




8 


10 


12 


500 


750 






100.0 


1.65(2) 


1.72(2) 


1.72(2) 


1.72(2) 


1.72(2) 
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FIG. 3: Interaction potential V^_^^{r, R,@) between H and II2 in a.u. The distance between 
hydrogen atoms in H2 is fixed at r = r(H2) = 1.409 a.u. 
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FIG. 5: Upper plot (a): total state-resolved cross section o"2-^o(v) vs. velocity v. Lower plot 
(b): corresponding thermal rate coefficients A;2_>o(T) vs. temperature T for the hydrogen molecule 
rotational relaxation process H2(j = 2 ^ j = 0) in H-H2 collision. 
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FIG. 6: Upper plot (a): total state-resolved cross section o"o^2(v) vs. velocity v. Lower plot 
(b): corresponding thermal rate coefficients A;o_>2(T) vs. temperature T for the hydrogen molecule 
rotational excitation process H2(j = — > j = 2) in H-H2 collision. 
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FIG. 7: Upper plot (a): total state-resolved cross section o"4_^2(v) vs. velocity v. Lower plot 
(b): corresponding thermal rate coefficients /c4_>2(T) vs. temperature T for the hydrogen molecule 
rotational relaxation process H2(j = 4 ^ j = 2) in H-H2 collision. 
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FIG. 8: Upper plot (a): total state-resolved cross section o"4_^o(v) vs. velocity v. Lower plot 
(b): corresponding thermal rate coefficients A;4_>o(T) vs. temperature T for the hydrogen molecule 
rotational relaxation process H2(j = 4 ^ j = 0) in H-H2 collision. 
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